
This file folder includes program for presenting trace plots of estimates in the MCMC sampling under semiparametric fitting with G=100 for cocaine use data. i.e. Fig.3 in Bayesian analysis of two-part nonlinear latent variable model: semiparametric method (statistical modeling). Among them,

  1. The routine 'main_SMI.c'is implemented to estimate unknown parameter  under semiparametric fitting with three staring values.

  2. The staring values are saved in files start1.dat, start2.dat and star3.dat. 

  3. The resulting estimates are saved in the fileEst_n.dat. 

  4. The Matlab routine 'converg.m' is used to plot the trace plots of estimates.

  5. The figures 'converg_exam.fig'and   'converg_exam.eps' are the trace plots of the estimates under semi parametric fitting with G=100 for cocaine use data.
